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Abstract 



Data on muon pair production obtained by the OPAL collaboration at centre 
^c of mass energies near the Z peak are analysed. Small angular mismatches between 

the directions of the two muons are used to assess the effects of initial state pho- 
ton radiation and initial-final-state radiation interference on the forward-backward 
asymmetry of muon pairs. The dependence of the asymmetry on the invariant 
mass of the pair is measured in a model-independent way. Effective vector and 
axial- vector couplings of the Z boson are determined and compared to the Standard 
Model expectations. 
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1 Introduction 

Many experiments have studied the forward-backward asymmetry of muon pairs produced 
in electron positron annihilation fl|-|3| , motivated by its sensitivity to interference between 
the axial- vector coupling of the Z and the vector couplings of the Z and the photon, the 
clean signature of the muons, and the lack of complications from t-channel exchange. The 
asymmetry depends on the centre of mass energy and, at tree level near the Z peak, this 
dependence is described by a straight line [|J, with a slope and intercept directly related 
to parameters of the Z boson. 

In the conventional method of analysis, however, the asymmetry is measured within 
a kinematic phase-space which integrates over the spectrum of radiated photons. This 
integrated asymmetry is then compared with the predictions of a theoretical model, whose 
parameters are varied to produce the optimum agreement. The integration over the 
spectrum of radiated photons noticeably changes the energy dependence of the asymmetry 
through two effects: 

• Initial State Radiation (ISR) lowers the effective centre of mass energy of the event, 
so the muon pair has an angular distribution which is appropriate to a lower energy 
and, furthermore, is distorted by the Lorentz boost. 

• Interference between photons emitted in the initial and final state (IFI) distorts the 
angular distribution from the usual 1 + cos 2 9 + a cos 9 dependence and produces 
a forward-backward asymmetry, strongly dependent on angular cuts, even in the 
absence of any axial coupling |J . 

These effects depend on the centre of mass (CM) energy (e.g., the intensity of ISR increases 
strongly when the energy exceeds the Z mass) and significantly distort the asymmetry at 
the levels of precision obtained by the LEP 1 experiments [fj. Q . 



In conventional analyses ISR is modelled by folding a radiator function calculated in 
QED with a Breit-Wigner model of the resonance cross section. For the effects of IFI, 
conventional analyses rely on the large cancellation expected if the cuts are set wide, i.e. 
when the cross section is integrated over a phase space which accepts almost all radiation 
from the initial and final states and the interference between the two |J. This integration 
represents a loss of angular information and the cancellation has not previously been 
verified experimentally. 

We present here an analysis which explicitly studies the effects of electromagnetic 
radiative corrections for muon pairs: 

• We consider ISR on an event-by-event basis. This enables us to assess the effective 
invariant mass of the muon pair, and thus study the variation of the asymmetry 
with energy, even at a single energy of the colliding beams. 

• We measure the effect of IFI on the angular distributions. By analysing small angu- 
lar mismatches between the directions of the two muons it is possible to identify and 
study those areas of the phase space where the IFI-induced asymmetry is significant, 
and see how the above-mentioned cancellation works. 

Our approach enables us to measure the Z boson effective couplings in a model- 
independent way, using only asymmetry measurements around the Z peak. The only 
assumptions used are those of QED, electron- muon universality, and the spin-1 nature of 
the Z boson. Our analysis is based on the measured properties of the [i + fi~ pair alone, 
and does not directly involve the detection of the radiated photons, which often escape 
direct detection either because of their low energy or because they go down the beam 
pipe. It thus includes the effect of soft radiative corrections, in contrast to approaches 
utilizing detected photons, which probe only hard corrections M. 



2 Choice of variables 

At tree level, muon pair production in e + e~-annihilation is a straightforward 2^2 
process: 

e + e~ — > fi + ~ \f~ (1) 

with the two final-state muons exactly back-to-back in the centre of mass (CM) system. 
For unpolarised beams, the azimuthal orientation of the event does not carry any useful 
information, and there is only one angle of interest, which can be chosen to be the polar 
angle of either of the two final-state muons 9 ± . Here we adopt the usual coordinate 
conventions, where the electron beam direction is along the +z axis of a right-handed 
cartesian system. The polar angles of the outgoing \i~ and /i + with respect to this 
direction are respectively 9~ and 8 + , and the corresponding azimuthal angles <p~ and 4> + 
are measured with respect to the x-axis, which points to the centre of the LEP ring. 

Higher-order corrections give rise to the radiation of real photons: 

e+e" -^ ^-7(7 . . .) (2) 



The radiated photons are not always directly observable: they may have very low energy, 
or may be radiated along the beam pipe, thus missing the detector altogether. But 
the photons can still be accounted for by measuring the angular mismatch between the 
directions of the muons. Since the two muons are no longer back-to-back, one needs three 
non-trivial angular variables to describe their directions. An especially convenient set is 
9* , r] and £. 

Following || the angle 9* is defined by: 

sing- - g+) 

cos ^ = -r-^ r-^r- (3) 

sin0- + sin0+ 

It reduces to cos#* = cos# + = — cos# - when 9 + = it — 9~ . In the case of ISR collinear 
with the beam direction, cos 9* equals cos 9 + in the CM system of the muon pair. 

The variable rj also depends only on 9 + and 9~: 

r, = |sin(g+ + ni (4) 

' \sm( y 9+ + 9-)\ + sm9+ + sm9-' v; 

If a photon is radiated exactly along the beam direction, r\ is the boost parameter and 
the energy of the photon, E^, is 

E-y = V y/s- (5) 

In this case, r\ measures the mismatch between the polar angles of the two muons, and is 
equal to zero if the muons are back-to-back. The invariant mass squared s' of the muon 
pair is then given by [] 

s'=(l-2r])s. (6) 

The acoplanarity £ is defined in the x — y plane: 



£ 



7T 



(7) 



It measures the angular mismatch between the two muons in the transverse direction, and 
is also equal to zero if the muons are back-to-back. 

Conventional analyses often use the acollinearity angle ( defined through the 3- 
momenta p± of the final muons: 

COsC = - p; rp; | - (8) 

It combines the angular mismatch in 9 with the mismatch in 0. This is not sensible in 
experiments where the resolution of the <fi measurement is significantly better than that 
of 9. Moreover, a mismatch in 9 is mostly due to strong ISR with the photon going along 
the beam direction, while a strong mismatch in is mostly due to photon radiation from 
one of the final state muons (FSR). These two processes, and their interference (IFI), 
have significantly different angular dependences, and their separation is an essential part 
of the present analysis. For these reasons, the acollinearity angle ( is not used here. 

^^This is true for any number of photons as long as they all have the same direction along the beam. 



3 Theoretical treatment 



3.1 Tree level formulae 



Consider first the (unphysical) case with no radiated photons, to which corrections will 
then be calculated. The normalised angular dependence of muon pair production at a 
CM energy yfs is given by: 



d 3 a(s) r 3 



| (i + C os 2 fl) + A FB (s) costf} 5(0 8{ri) (9) 



a(s) dcosO d£ dr] I 8 

which has a trivial dependence on £ and 77, described by the two Dirac ^-functions. The 
coefficient Afb(s) of the term linear in cos 9 is the forward-backward asymmetry: 

/ \ <j(cos9 > 0) — cr(cos9 < 0) , . 

Afb{s) = a(cos9>0) + a(cos9<0) (10) 

= **« (11) 

Assuming only that the Z boson is a massive spin-1 resonance with mass M z and width 
Tz, and following the notation of ||, the energy- dependent functions ^1,3(5) near the Z 
peak have the following form [|3],|9|]: 

F 1 (s) = l + 2Re{x*o(s)C; z }+\xo(s)\ 2 C s zz, (12) 

F 3 (s) = 2Re{ X * (s)C» z } + 4\ X o(s)\ 2 C a zz , (13) 

where 

MS) = ICs-Ml + iMzTz (14) 

Assuming electron-muon universality in vector and axial-vector couplings, gy A = gy A = 
gv,Ai the coefficients C take the form: 

C°z = 9v , C zz =(g v + g A )\ (15) 

C a lZ = 9 2 A , C zz = g 2 v g 2 A- (16) 

In the immediate vicinity of the Z pole the asymmetry (|TT|) is a linear function of s: 

* , s ZC a zz s-M z T ^ 3Ci. z 

A fb(s) = —p. + ^^K—^ (17) 

U ZZ z " b ^zz 

3g 2 v g\ s - Ml v 3g\ 



(9 2 v + 9 2 a) 2 + 2, K (g 2 v + g 2 A r (18) 

The first (constant) term in both ([17]) and (|18"), the asymmetry at peak, depends only on 
the ratio gv/gA, while the slope with energy of the second term allows one to measure the 
axial-vector leptonic coupling, g A . The constant K, which stands for the ratio of the Z 
boson and photon propagator normalisation factors, determines the scale of the gy and g^ 



parameters: only two of the three quantities gy, gA and K are independent, if asymmetry 
is the only measured quantity. 

Note that eqs. fllTD and (|18D do not contain Y Z - In particular, the slope of the asym- 
metry with energy is independent of V z- The linear approximation remains valid over the 
region where the Z dominates the photon. 

In the Standard Model, the constant K is expressed through the ratio of the electro- 
magnetic coupling a and the Fermi constant Gp\ 



G F Mj 



'< - %?£ (w) 



After the imaginary part of the photon propagator is taken into account, which results 
in a small offset in the forward-backward asymmetry at peak [12|], eqs. ( |12| — [19]) still 



hold [pTOj pT2[| in what is called the "Improved Born Approximation" , with gy and g_& now 
standing for the real parts of the effective vector and axial- vector couplings of the Z boson. 
The contributions of the imaginary remnants of the effective couplings are small, and 
have been neglected in this analysis. The numerical value of the normalisation constant 
K changes as the scale dependence of the electromagnetic coupling is taken into account: 

a -► a(M|) « 1/128.89. (20) 

3.2 Radiative QED corrections 

Initial state radiation (ISR), final state radiation (FSR) and the interference of the two 
(IFI) affect the angular distribution of final muons in different ways. 

ISR photons are radiated mainly along the beam axis, and the CM frame of the muon 
pair acquires a boost. As mentioned above, under collinear ISR 9' remains equal to 9 + 
in this frame. The acoplanarity variable £, defined in the transverse plane, is also largely 
insensitive to ISR. The parameter 77 however, is essentially proportional to the energy of 
the emitted photon. Thus, events with significant ISR typically have 77 significantly larger 
than £. The angular distribution 



d 2 a 



dcos9' (h) 



{ ^ (l + cos 2 6')+A FB (s') cose'}f(r]) (21) 



acquires a non-trivial ^-dependence described by the function f(r]), and an additional 
^-dependence through the argument s' — (1 — 2rf)s of the asymmetric term Afb(s'). So, 
the measurement of the mismatch in the polar angles of the two muons, described by 
the variable 77, allows the forward-backward asymmetry A F b to be measured directly at 
various energies VV, below and up to the actual initial CM energy, y/s. 

FSR is essentially symmetric around the final muon direction, and an angular mis- 
match in the longitudinal direction is close to that in the transverse direction, yielding on 
average 77 ~ £. FSR is mainly directed along the final muons, and its effect alone on the 
forward-backward asymmetry is unmeasurably small. 

Radiation with significant initial-final interference, IFI, is concentrated also mainly in 
the areas where r) c± £. It is a complicated function of all three angular variables and 



contains a term which is odd in cos#*, thus introducing an additional forward-backward 
asymmetry. This is expected to be positive for softer photons, E 1 < Tz/2, and negative 
for harder ones. This strong variation of the asymmetry as a function of acoplanarity 
has been observed in jn| . For a totally inclusive treatment (when one integrates over all 



values of r) and £), IFI should have a negligible effect on the asymmetry, as these two 
regions cancel each other almost completely ||. By measuring t] and £ we can identify 
those areas of the phase space where this IFI-induced asymmetry is significant and check 
the mechanics of the expected cancellation. 

Consider a tight cut on acoplanarity £ < £ - This restricts the phase space of the 
radiated photon and disturbs the delicate balance required for the cancellation of the 
initial-final interference effects. The additional asymmetric term, which enters the ex- 
pression for the full differential cross section, has a characteristic logarithmic dependence 
on cos 0' 00: 



f* J^rhr, ~ i I (1 + cosV) + AfM1 ~ 2 " )s) cosS '} 

( 1 1 _i_ pns ft* i 

x{/ + (^)[l + ^,cosg')] + /-(^ln 1 _^ g . }. (22) 

The functions (3(r), cos 9') and f±(rj) implicitly depend upon the value of the acoplanarity 
cut £o! a complete theoretical calculation of these functions for arbitrary values of the 
cut £o would involve a detailed analysis of the participating interfering diagrams and the 
complete set of the SM parameters. However, the treatment is significantly simplified if 
the cut £o is chosen to be very tight, of order 10~ 3 . In this case, the functions f± and (3 
can be calculated by numerical integration of the analytic formulae describing QED ra- 
diative corrections in the single soft photon approximation [TJ[]. By comparing theoretical 
formulae and generator-level Monte Carlo simulations with [|Tj| 16|] and without [1_1. 17 



exponentiation of radiative corrections, it was verified that the above approximation is 
adequate at energies around the Z peak, if £ and rj are small enough and one does not get 
too close to the edges of cos 6* distribution. 

In this approximation, the term with the function f+{rj) contains the contributions 
from ISR, which is cos#'-independent, and from the remnants of FSR, which result in 
an even cos ^'-dependence, described by the function /3(j], cos#*). The latter can be 
parametrised as 

R( V\ f-W 1 1 1+COsg * fn* 

pin, cos ) = ——— - — In- — (23) 

w ' /+(?/) 4cos#* l-cosfl' 

and, in the range of the angular variables considered in this analysis, represents a small 
(~ few percent) correction to the main symmetric term. 

The last term in (P2"D, with the function f-{f]), describes the contribution of initial-final 
radiation interference, and has an odd cos ^'-dependence. This term is responsible for the 
additional asymmetry which arises because of the tight acoplanarity cut. Fortunately, 
the cos ^'-dependences of the two odd terms in fl2"2"|) are significantly different, and, given 
the statistical power of the LEP 1 data set, these two asymmetric contributions can be 
separated. 

8 



4 Monte Carlo samples 

Three Monte Carlo samples have been used at various stages of this study. They were 



obtained using KORALZ generator version 4.0 |L7| and full OPAL detector simulation 
It is an important feature of this analysis that the numerical results obtained in Section 
[7] are, in fact, almost independent of these simulations, the samples being used only to 
assess systematic errors and derive a few small corrections. 

The first sample contains 600000 muon pair events with soft photon exponentiation, 
but without the initial-final radiation interference. This sample, labelled "MCI", was 
used to study resolutions and efficiencies, as described in subsections 5.2 and 7.1. 

The second sample, labelled as "MC2", contains 100000 events with the complete set 
of 0{ot) corrections, which includes the initial-final radiation interference. Alongside with 
MCI, MC2 is used in Section 6 to illustrate the differences in the angular distributions of 
the various data subsamples. 

The third sample contains 600000 r pair events. It is referred to as "MC3", and is 
used for background studies together with MCI, as described in the following section. 

5 Event selection 

Events collected by the OPAL experiment at LEP 1, at and around the Z peak during 
1993, 1994 and 1995 are used in this study. The data correspond to total integrated 
luminosities of about 82 pb _1 at the CM energy around 91.22 GeV (hereafter referred to 
as "peak"), 17 pb^ 1 at 89.44 GeV ("peak-2") and 18 pb' 1 at 92.97 GeV ("peak+2"). 

A detailed description of the OPAL detector is given elsewhere [H],|2(J . Here we briefly 



describe some subdetectors relevant to our analysis. Most of these form a set of coaxial 
cylinders, with varying coverage in 9. The innermost is the central vertex detector, which 
consists of a set of axial (CVA) and stereo (CVS) wires, both with angular coverage 
| cos#| = 0.92. This is followed by the multi-wire jet chamber (CJ), which is surrounded 
by thin ^-chambers (CZ), the latter covering the angular range | cos 6*| < 0.72. The 
coaxial magnetic field of 0.435 T allows the measurement of the momenta of charged 
particles. The following layer is formed by the electromagnetic (ECAL) and hadronic 
(HCAL) calorimeters. The outer layer consists of a set of barrel (MB) and endcap (ME) 
muon chambers. 



5.1 Selection of Muon Pair Events and Background Rejection 

Muons are identified either by hits in muon chambers (MB and/or ME) or by the specific 
pattern of energy deposition in HCAL or ECAL. Events passing the multihadron and 
cosmic veto and containing 2 tracks identified as muons are selected if the visible energy 
in the event (defined as the sum of the two muon energies and the highest energy ECAL 
cluster in the event) is larger than 0.6-y/i (see [|[ for details). This part of the e + e~ — > 
fi + fi~ selection is identical to that described in ||, except that here we do not require the 
cut on the acollinearity angle (||), but do require that the two muon tracks be measured 
to have opposite charges. 



The visible energy cut is designed to reduce the only significant background to the 
processes (|1|,@), t pair production, studied using Monte Carlo samples MCI and MC3. 
In the usual inclusive asymmetry analysis this cut removes most of the tt background, 
reducing its fraction in the data event sample to about 1%. However, tt events have a 
much broader 77-distribution than \x\x events, and even this small amount of background 
can interfere with the analysis presented here. An additional cut is applied to the missing 
transverse momentum in the event: 



x T = y ( s P*) 2 + (^) 2 < o.l, (24) 

where the sums are taken over respective components of the momenta of the two charged 
particles and the two most energetic electromagnetic clusters in the event. The cut ([24]) 
removes ~ 75% of the remaining background, as shown in fig. [I]. For the asymmetry 
analysis described in Section 7, we also impose a strong requirement on the acoplanarity, 
£< 0.004, restrict | cos^'l < 0.92, and limit s' to the region where the asymmetry remains 
the linear function of energy by imposing rj < 0.06,0.08 and 0.10 at the peak— 2, peak 
and peak+2 energy points, respectively. 

These additional cuts eliminate the charge- dependent tracking problems described 
in || , which lead to asymmetry biases in the subsample of poorly measured events retained 
for the conventional OPAL asymmetry measurements, and further reduce the r pair 
background. The ^-distribution is much broader for the r pair background than for the 
\i pair signal, so the tight cut on acoplanarity rejects most of the remaining background 
events. For the asymmetry sample, used in our fits in Section 7, the estimated number of 
remaining tt events is 25 ± 3 out of 66143 selected events. 

As described in the next subsection, we also require that all tracks in the asymmetry 
sample be measured by subdetectors with optimum resolution, which retains about 62% 
of otherwise selected events. 



5.2 Angular resolutions and muon pair event classification 

It is essential for this analysis to achieve the best possible angular resolutions in both 
9 and <fi. For this reason, both muon tracks are required to have at least one hit in 
the central vertex detector in both axial (CVA) and stereo (CVS) planes, giving a good 
measurement of the production vertex. This requirement rejects about 25% of selected 
events. The z-resolution of the central jet chamber (CJ) is insufficient for this analysis, and 
therefore measurements are required from either the z-chambers or the muon chambers. 
This requirement rejects a further 13% of the selected events. 

The azimuthal angle measurement is made by the jet chamber CJ, which determines 
the resolution in £. A graphical illustration of the effect of the finite £ resolution is 
given in fig. ^|a, which shows the distribution of a subsample of events from the central 
region, | cos^*! < 0.2, with respect to the variable log(l/£). With a perfect detector, this 
distribution is expected to fall slowly and smoothly after reaching a plateau at log(l/£) ~ 
2.5 — 3. With a finite resolution og, all events with unmeasurably small true angular 
mismatch £ <C a^, which would have appeared at very high values of log(l/£) if the 
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measurements were perfectly accurate, acquire larger £ of order of a^, and gather into the 
sharp peak at log(l/£) ~ log(l/cr^) ~ 3.5. The solid line in fig. ||a represents the result 
of a Gaussian fit (in £) to the data points, resulting in the estimated resolution for this 
subsample, a^ = 0.36 mrad, with a statistical error about 1%. The £ resolution remains 
virtually constant in the barrel region of the detector, |cos#*| < 0.72, and increases 
gradually to ~ 0.85 mrad for |cos#*| > 0.85. 

The resolution of the r\ measurements depends on the sub detector used to measure 
the far ends of the two tracks in the event. CZ provides the best measurement, followed 
by MB, ME and, finally, CJ. The selected muon pair events are classified according to the 
effective resolution of r\ measurement: 

1. Polar angles of both tracks are determined from hits in the CZ. These events have 
the best resolution in 9. 

2. Muon chambers are used to measure both polar angles. 

3. One of the tracks has its 9 determined from the CZ measurement while the other is 
measured with the muon chambers. 

4. Either or both tracks have their polar angle determined by CJ, and/or have no hits 
inCVA/CVS. 

The distribution in log(l/?7) of a subsample of class 1 events from the same region 
| cos#*| < 0.2 is shown in fig. |2|b. Here too, the position of the peak defines the resolution 
of the detector. For this particular subsample, the Gaussian fit (in rj), shown by the curve 
in fig. ||b, gives a v = 0.95- 10~ 3 , with a statistical error of about 1%. Like the £ resolution, 
the f] resolution increases for cos 9* values outside the barrel region. 

Both £ and 77 resolutions were determined separately for 10 bins of |cos#*|. Non- 
Gaussian tails were approximated by a Breit-Wigner distribution, with parameters deter- 
mined from the Monte Carlo sample MCI. The r\ resolution was determined separately 
for the different event classes defined above. Numbers of selected events in each class are 
shown in table [I] together with respective estimated resolutions. 

Only events from classes [l] and || are used in our main analysis (adding up to about 
62 % of the selected muon pair events remaining after the cut on missing transverse 
energy), with class |3| events used for systematic studies. The r] resolution for events from 
class [| was found to be non-uniform with respect to the polar angle and generally poor, 
so that they could not be used. 



6 cos 9* distributions for various 77, £ regions 

The scatter plot of the class [1] events in the log(l/?y) — log(l/£) plane is shown in fig. [3[ 
As explained above, the events placed along the diagonal of this plot, corresponding to 
T) ~ £ (area with dense horizontal shading), have a high probability of strong FSR and/or 
IFI. Events with rj ^> £, scattered above the diagonal (diagonally shaded triangular area 
in fig. [5p typically have a high probability of significant ISR. As in fig. |2|, the position of 
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(E cm ) 


Class 1 


Class 2 


Class 3 


Class 4 


Total 


p-2 


89.45 GeV 


3215 


1265 


579 


2082 


7141 


peak 


91.22 GeV 


46727 


18858 


9401 


31244 


106230 


p+2 


92.97 GeV 


4853 


1896 


889 


3147 


10785 


Total 




54795 


22019 


10869 


36473 


124156 




Resolutions 














(T V , 10- 3 


0.95 


2.0 — 3.4 


1.5 


2.4-- 8.0 






<j£, mrad 


0.35 


0.35 — 0.85 


0.35 


0.35 — 2.5 





Table 1: Numbers of events (prior to any cuts on angular variables cos6', rj and £, but 
after all other cuts) and estimated resolutions in rj and £ for the three energy points and 
various classes of events as defined in the text. Events from classes 1 and 3 are from the 
barrel region of the detector, | cos^'l < 0.72; most of class 2 events belong to the region 
0.72 < | cos^'l < 0.92. For classes 2 and 4 the resolutions depend on the polar angle: the 
smaller numbers for the resolutions a v , a^ correspond to | cos^*| > 0.72, while the larger 
numbers refer to the edge of acceptance | cos^'l < 0.92. 



the peak in fig. |3] is determined by the detector resolution in rj and £. Class 2 events show 
a very similar scatter plot, apart from the fact that, because of the inferior resolution, the 
peak is shifted towards higher values of £ and rj (i.e. down along the diagonal of the plot 
in fig. |). 

A cut on acoplanarity £ = 0.004, corresponding to log 10 (l/£) ~ 2.4, is shown in 
fig. by the horizontal dashed line. The vertical dashed line corresponds to rj = 0.008, 
log 10 (l/?7) ~ 2.1. This value was chosen to separate the area of relatively large rj values, 
where the finite ^-resolution effects are not too important, and the region of small r), 
where the distributions are significantly smeared by the detector resolution. 

The upper-left quadrant in fig. ^ is filled with events with small £ (£ < 0.004) and 
large rj [r\ > 0.008). In this kinematic range there is a high probability of strong ISR, but 
it is essentially free of IFI and FSR contributions. Thus the distribution of the events 
from this quadrant in cos 6* should be well described by the usual quadratic function of 
cos#', eq. ([21]). 

The upper right quadrant in fig. |3] contains events with both £ and rj small. Here the 
IFI contribution is expected to be strong and should lead to an additional logarithmic 
dependence on cos#", described by the second term in curly brackets in eq. (22). This 



dependence gives rise to an additional, IFI-induced forward-backward asymmetry, which 
is expected to be positive in this part of the r)-£ plane. 

The lower left quadrant of fig. § also contains a strong IFI contribution, but in this 
region the IFI-induced asymmetry is expected to be negative. Theory predicts that, when 
integrated over the whole range of £ and rj, these positive and negative logarithmic inter- 
ference terms cancel each other so that the overall cos#* distribution is again described 
well by a simple quadratic function of cos 6', with the only asymmetric term being linear 
in cos 6*. 
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Figs. [||a-d show the cos ^'-dependence for class 1 and class 2 events collected at the Z 
peak. The measured distributions were corrected for efficiency and background using the 
Hjj, and tt Monte Carlo samples including full detector simulation. Figs. £|a-c correspond 
to upper left, upper right and lower left quadrants of fig. [3] respectively, while fig. f|d shows 
the cosR' distribution integrated over the whole 77— £ plane. Also shown are fits to the 
measured distributions using the following function: 



1 1 4- ms ft* 

{a + b cos R' + c cos 2 ft') (1 + d- In -— n ). (25) 

2 1 — cos R* 

which is merely a simplified version of eq. fl2"2"|) after the integration over the respective 
quadrant of the 7/-£ plane. Among the four fit parameters in (p5j), a,b,c and d, the 
last one, d, is particularly interesting, as it determines the amount of the IFI-induced 
forward-backward asymmetry. Its fitted values are shown in the figure for each of the 
four distributions. 

These fit results clearly demonstrate that theoretical expectations are fulfilled: the 
cos 6**-distribution of the events taken from the upper left quadrant of fig. |] is indeed well 
described by a quadratic function, giving d = 0.01 ± 0.07, compatible with the absence 
of IFI contribution. On the contrary, events from the two quadrants situated along the 
diagonal of fig. || both show a significant non-zero IFI contribution: the upper right 
quadrant yields a positive IFI-induced asymmetric term, d = 0.09 ± 0.03, while the lower 
left quadrant (events with both £ and rj large) yields a negative value, d = —0.18 ± 0.05. 
Most notably, when all regions of T) and £ are summed, the resulting cos R* distribution 
in fig. [p, with d = 0.03 ± 0.03, is again compatible with the absence of the IFI-induced 
logarithmic term. 

Before moving on to the detailed quantitative analysis of the measured angular distri- 
butions, let us show that the differences in the asymmetric parts of the measured cos R' 
distributions in different r)—£ regions are indeed caused by the initial-final radiation inter- 
ference. Consider the "differential forward-backward asymmetry", defined by 

dcr(cos8 m ) da{— cos 6*) 

A frmcO'\ — dcosO* dcos9 m [1d\ 

A FB (cosft ) = dff(cose .) dff( _ -7^T > ( 26 ) 

dcos8' ' dcose* 

as a function of cos6**, for the same three separate areas of the log(l/£) — log(l/?y) plane 
considered above. This is compared to two Monte Carlo samples, MCI and MC2, as 
defined in Section 4. MCI has the IFI term switched off, while MC2 includes the IFI 
contribution. 

Fig. [5]a represents the upper left corner of fig. ^ rich with ISR, and shows a negative 
asymmetry with the smooth angular dependence ~ cos#*/(l + cos 2 ^*), typical of the 
linear term. The two lines, corresponding to the Monte Carlo samples with and without 
initial-final interference, display no significant differences, and both describe the data well. 

In contrast, the shape of the angular dependence for the upper right (fig. |5|b) and lower 
left (fig. |5|c) quadrants is dominated by the logarithmic term ~ ln[(l+cos0*)/(l— cos0*)], 
characteristic of the IFI-induced forward-backward asymmetry. The data and the MC2 
sample agree reasonably well, while for the MCI this is not the case. 
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When integrated over the whole r/-£ plane, the angular dependence is essentially sym- 
metric, as shown in fig. |5|d, without any significant deviation between the data and either 
Monte Carlo sample. 



7 Asymmetry analysis 

7.1 Probability density 

As mentioned above, a tight cut on acoplanarity selects an area of phase space where 
we can apply the formalism presented in section ^|, measure the energy dependence of 
the forward-backward asymmetry and determine the vector and axial- vector couplings of 
Z. Based on the expression (^2|) for the double-differential cross section, one obtains the 
following probability density function: 



V(t), cos 9') ~ { - (1 + cos 2 6') + 
I 8 



A + A' £ 



cos 9' 



xW ,){ 1+%ros , )+B |M> h i±^;} (27) 

where Aq, A' and B are constants to be determined from a fit to the data. Aq and A' 
correspond to the coupling combination terms in equation (|18|). In the single soft photon 
approximation of QED one expects B — 1, so by measuring B we can accommodate and 
measure deviations from this approximation. 

The function f+(rj) essentially defines the shape of the ^-dependence of the cross sec- 
tion, and was measured directly from the data. Indeed, the single-differential distribution 
with respect to 77, obtained by integrating the probability density ( P7| ) over the whole 
cos#* range, is essentially proportional to the function /+(?/). The corrections from the 
FSR contribution f3, defined in eq. (E3j), and the product of the two terms in eq. (P7|) 
which are odd functions of cos6' are fairly small and can be easily taken into account. 

Before application to the data, rj resolution smearing must be explicitly applied to 
(P7p. Note that the rj dependence enters not only through f±, but also through s' within 
the first pair of curly brackets, so the probability distribution folded with the resolution is 
no longer factorisable. However, the //-dependence of the first term is linear, so one only 
needs to calculate three different ratios of folded functions: ^i(r/) = /_(r/)// + (r/), ^(t?) = 
Vf-(v)/ f+(v) an d ^(t?) = vf+iv)/ f+(v)i where / denotes a function / convoluted with 
the ^-resolution. These three ratios are presented in fig. |6], together with similar ratios 
without resolution smearing. Even the largest of these three ratios, \l/i = /_//+, which 
essentially determines the IFI contribution to the asymmetry, is much smaller than 1 at 
rj ~ and quickly becomes even smaller outside a narrow range of rj, the width of which 
is governed by the cut on £ and the resolution in rj. 

Since the standard OPAL selection efficiency for /i-pair events is very close to 100%, 
the efficiency of the class 1 and 2 selections can be determined directly from the data. 
Fig. shows the ratios, £j(cos0 # ), of class 1 and class 2 events to events of all classes, as 
a function of cos#*, for events passing all other requirements for the asymmetry sample, 
including cuts on xt, £, V and cos#V The small inefficiencies due to the xt requirement 

14 



and resolution losses in the £ cut are calculated using the Monte Carlo, as described in 
the Appendix. 

Finally, the expected probability distributions (^7p summed over event classes, con- 
voluted with the respective r\ resolutions and appropriately weighted with the respective 
efficiencies, are normalised so that the total probability is independent of the fit parame- 
ters. The resulting formula, incorporating all the corrections described above, is given in 
the Appendix. 

It is convenient to re-express the asymmetry at peak, Aq, and the slope of the asym- 
metry with energy, A', in terms of the vector and axial-vector couplings of the Z, as in 



18) , and use the following set of fit parameters: 



Pi= l H P 2 = %=, V, = B. (28) 



\9a\ " \/K 

The first two parameters have obvious physical meanings, while the last determines the 
measured intensity of the IFI-induced term, compared to the single soft photon approxi- 
mation. The asymmetry at peak and the slope now read: 

A = A FB (M 2 z )=3p 2 1 /(l+p 2 1 ) 2 + A^ (29) 

A' = AUMf)=3(l/p 2 ) 2 /(l+Pi) 2 , (30) 

where A^ = 0.002 is the offset to the pole asymmetry due to the imaginary part of the 
photon propagator. 

7.2 Results 

An unbinned maximum likelihood fit is made using the data sample from classes [l] and 0, 
with the probability density function and the set of parameters described in the previous 
subsection. The cut on acoplanarity is chosen to be £o — 0.004, which is found to be small 
enough to reject most of the FSR contribution and justify the soft photon approximation, 
while simultaneously being large enough compared to the £ resolution. The range of cos 9* 
used in the fit is limited by the acceptances of relevant subdetectors to —0.92 <cos9'< 
0.92, while the upper bound for the variable r], r] majX , is limited by the range of s' where the 
asymmetry is expected to be a linear function of energy. We choose r] max = 0.06, 0.08, 0.10 
for data taken at peak— 2, peak and peak+2, respectively. This upper bound for rj removes 
about 0.1% of the remaining events. 

Table presents fit results for peak data only, for peak— 2 and peak+2 simultaneously, 
and for all three energy points simultaneously. The errors shown in the table are statistical 
only. The corresponding correlation matrices are given in table 0. No meaningful results 
have been obtained for peak— 2 or peak+2 data sets separately, because of the limited 
statistics at these energies. 

The unbinned maximum likelihood fit does not give any goodness-of-fit parameter for 
judging the quality of the fit. In order to do this and to illustrate our results graphically, 
we subdivide the data into 30 bins in v?, and perform a single parameter maximum 
likelihood fit in each bin (with B fixed to its previously determined value, B = 0.840) for 
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E 


\9v/9a\ 


\9a\/VK 


B 


Z peak 


0.0830 ±0.0125 


0.6013 ±0.0314 


0.811 ±0.129 


p-2 & p+2 


0.0720 ±0.0312 


0.6226 ±0.0212 


0.952 ±0.322 


All energies 


0.0795 ±0.0114 


0.6165 ±0.0177 


0.840 ±0.120 



Table 2: Fit results for various energy combinations. Errors are statistical only. 



the coefficient of the cos#* term in eq. (ffi\ )- The results are presented in fig. |8|. One sees 
that the measured asymmetry at various v? values are indeed aligned close to a straight 
line, whose value at s' = M§ and slope can now be determined from a minimum \ 2 fit 
to these points. The fitted line is also shown in fig. ||. The value of x 2 /d.o.f.= 38.9/28 
suggests that the fit quality is acceptable. The results of this fit: 



9v 

9A 



0.0813 ±0.0082, 



9a 



0.6246 ±0.0184 



(31) 



are in agreement with the results of our main fit from table [|. The smaller error in (|31|) 
is due to the fact that the parameter B was fixed; fixing B in the maximum likelihood fit 
also results in smaller errors, 0.0082 and 0.0175, respectively. 

7.3 Systematic Studies 

The probability density function used in the fit depends upon a number of parameters 
whose values cannot be precisely fixed. The variation of fit results due to varying these 
parameters within a reasonable range allows one to estimate corresponding systematic 
errors. 

Various sources of systematic error have been considered, and the resulting errors are 
summarised in table |], for the data taken at the Z peak only, and for all three energies 
analysed simultaneously. 

1. The fit was repeated with an additional cut Icos^'l < 0.90, to check sensitivity 
against the variation of the edge of the geometric acceptance. The assigned sys- 
tematic error is the absolute value of the shift, wherever the shift is statistically 
significant, plus a small contribution due to the uncertainty of the absolute scale of 
the cos^* measurement. 
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Z peak 


\9v/9a\ 


\9a\/VK 


B 


\9v/9a\ 


1.000 


-0.396 


-0.698 


\9a\/VK 


-0.396 


1.000 


0.190 


B 


-0.698 


0.190 


1.000 


p-2 & p+2 


\9v/9a\ 


\9a\/VK 


B 


\9v/9a\ 


1.000 


-0.186 


-0.663 


\9a\/VK 


-0.186 


1.000 


0.113 


B 


-0.663 


0.113 


1.000 


All energies 


\9v/9a\ 


\9a\/VK 


B 


\9v/9a\ 


1.000 


-0.261 


-0.694 


\9a\NK 


-0.261 


1.000 


0.133 


B 


-0.694 


0.133 


1.000 



Table 3: Correlation matrices for the three fits. 



2. The parameters a v , describing the experimental resolution in the r\ measurement, 
and determined from the data in bins of cos 9* for various classes of events, were 
scaled by a factor of 1 ± 0.1 for each class separately. The assigned systematic error 
is the largest of the absolute values of the observed shifts. 

3. The calculations involving the function f-(rj) are less reliable for 77 > Tz/(2Mz), 
where the photon spectrum can be affected by the Z resonance lineshape. To study 
the influence of this uncertainty, the fit was repeated with f-(rj) set to zero for 
77 > 0.010. The assigned systematic error is the absolute value of the shift. 

4. In order to check for possible biases due to the approximations made in deriving 
eq. (|i~8|) , the fit was repeated with next-to-leading terms in 77 taken into account. 
The assigned systematic error is the absolute value of the shift. 

5. Monte Carlo studies have shown that deviations from the equation s' = s(l — 2rj) 
within the angular range considered here do not exceed ±0.5%. Possible biases 
were checked by replacing 77 with 8\ + (1 + #2)77, where 5\^ = ±0.5%. The assigned 
systematic error is the largest of the absolute values of the shifts. 

The total systematic uncertainty for each fit parameter was calculated as a quadratic 
sum of the partial contributions. 

The following checks have also been made: 

• The fit was repeated with an additional cut to remove data in the range 0.70 < 
I cos #*| < 0.75, to exclude the edge of the barrel part of the detector. 

• The acoplanarity cut £0 was varied by ±0.001 from its central value of 0.004. 

• The upper limit of the 77 range was varied by ±0.010 from its central value for each 
energy point. 
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Variation 


Z peak 


All energies 


\9v/9a\ 


\9a\/VK 


B 


\9v/9a\ 


\9a\/VK 


B 


1 


| cos 6'\ < 0.90 


0.0014 


0.0000 


0.054 


0.0002 


0.0000 


0.062 


2 


a n 


0.0005 


0.0009 


0.002 


0.0006 


0.0002 


0.002 


3 


f-ij]) tail 


0.0001 


0.0038 


0.001 


0.0006 


0.0012 


0.003 


4 


s'-dependence 


0.0005 


0.0053 


0.004 


0.0007 


0.0018 


0.001 


5 


s' *->■ 7] relation 


0.0001 


0.0015 


0.001 


0.0001 


0.0016 


0.001 




Total syst. 


0.0016 


0.0068 


0.054 


0.0011 


0.0027 


0.063 




Stat, error 


0.0125 


0.0314 


0.129 


0.0114 


0.0177 


0.120 




Total error 


0.0126 


0.0321 


0.140 


0.0115 


0.0179 


0.136 



Table 4: Various contributions to the systematic errors, for the data taken at the Z peak 
only, and for all three energy points. 



The number of bins in the measured ^-dependence was changed by ±10 from its 
default value of 50. 

In order to check the reliability of the efficiency calculation, class 3 events (defined 
in subsection 5.2) were added to the analysis. 

The number of bins in the measured efficiency as a function of cos 9* was changed 
by ±50 from its default value of 100. 



• The cut on the missing transverse energy ( f24|) was tightened from 0.10 to 0.05, 
effectively reducing the r pair background contribution by a factor of 2. 

In all these cases, the observed shifts were well within expected statistical variations, 
each of which constituted a fraction of the total statistical error, therefore no additional 
systematic errors were assigned. 

7.4 Determination of gy, Qa and sin 2 9w 

Combining statistical and systematic errors, for the results at all three energy points we 
obtain: 



9v 



9A 

9a 



K 
B 



0.0795 ±0.0115 

0.6165 ±0.0179 
0.840 ±0.136. 



(32) 

(33) 
(34) 



These results are essentially independent of SM assumptions and SM parameter values. 
The only assumptions used are those of QED, electron- muon universality and the spin-1 
nature of Z. 
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For comparison, the values for the quantities fl3"2| ) and ([^J), extracted from the ratios 
C zz /C zz an d C^ z /C zz , as measured in the conventional asymmetry analysis || using 
the full OPAL muon data sample of 1990-1995, are: 



9v 



9A 

9a 



0.0713 ±0.0055 
0.6178 ±0.0147. 



(35) 
(36) 



These two sets of numbers are in agreement with each other, as well as with the SM 
expectations Q: 



9v 



9a 
9a 



n729+ - 0015 



0.59459; 



-0.00047 
-0.00013- 



(37) 
(38) 



The measured value of B, eq. (0), is also compatible with the expectation, B = 1, of 
the single soft photon approximation. 

From the measured ratio \gv/ 9a\ we can directly determine the effective weak mixing 
angle in the Standard Model (assuming that gy and g& have the same signs): 



sin 2 C = 



9v 

gA 



0.2301 ±0.0029, 



(39) 



which is in agreement with the world average 0.23150 ±0.00016 [|§. In order to determine 
the effective couplings gy and g& separately, we have to substitute numerical values (which 
are well measured elsewhere || in the context of the Standard Model) into the definition 
of the normalization constant K (eq. (|19|)). Using (p0|), one gets yK = 0.843108, which 
gives the following values for the vector and axial-vector couplings of the Z boson: 



\g v \ = 0.0413 ±0.0060, 
\g A \ = 0.520 ±0.015. 



(40) 



The higher precision of the conventional analysis is mostly due to its higher statistics, 
since this analysis is restricted to events with accurate angular measurements. However, 
in contrast with the conventional analysis, this analysis has the ability of extracting the 
slope of the energy dependence of the asymmetry (and hence the parameter g&) from the 
data taken at a single energy point. By comparing the errors on the parameter g^/yK 
determined from the peak and off-peak data in Table ^| one can see that the weight 
of the peak contribution to the final precision is quite significant. This information is 
clearly complementary to the standard analysis, and can be combined with the results 
of the latter to improve the overall precision on the relevant coefficient, C%. The values 
for these coefficients from the standard OPAL analysis ||, from muon data only, and 
averaged over the three lepton flavours, 



c^V 



0.232 ±0.011, 



C* z (l + r) = 0.2350 ± 0.0080, 



(41) 



.,(5) 



2 Assuming M H = 150±|g° GeV, M t = 175 ±5 GeV, a s = 0.119 ±0.002 and Aa { h J > = 0.02804 ±0.0065. 
See [II for details. 
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can be combined with the value obtained, using eq. ([H]), from this analysis (Z peak only): 

C^Gu+/0 = 0.257 ± 0.027. (42) 

We obtain: 

C^z^V - ) = 0.236 ± 0.010, C* z (l + r) = 0.2368 ± 0.0077, (43) 

which now represent the best OPAL values for these coefficients. 

8 Conclusion and outlook 

We have analysed the angular dependence of muon pair production in electron-positron 
annihilation at centre of mass energies near the Z peak, using various angular variables. 
Our approach is novel in a number of respects. The usual procedure involves the inte- 
gration over the phase space of the radiated photons, limited by a cut on acollinearity 
(eq. (|8p). In contrast, we measure small angular mismatches between the directions of 
the two final muons, separately in the polar (if) and azimuthal (£) directions, and use 
them to determine the influence of the initial and final state photon radiation and their 
interference. Effects of final state photon radiation are removed by applying a tight cut 
on the acoplanarity, £. The contribution of the additional asymmetric term arising as a 
result of this cut is measured through its specific polar angle dependence. The variable 77 
is used to assess the energy of the radiated photon and to determine the variation of the 
forward-backward asymmetry with the invariant mass of the muon pair, which is shown 
to be linear in the vicinity of the Z peak (see fig. ||). 

By using a well-behaved variable, cos 6', instead of the polar angle of one of the muons, 
and explicitly incorporating the initial-final interference into the fit, we significantly reduce 
the dependence of the measured asymmetry upon the polar angle acceptance cut. 

The measured values presented in equations (p2}-p4|) are directly obtained from the 
data; they can be compared to those of other experiments, or to theoretical models. 
By substituting the SM value for the constant K, we get results for g v , g& and the 
effective weak mixing angle compatible with those obtained with the analyses based on 
the assumptions of the Standard Model. The statistical precision of our result, while 
obviously inferior to that of the model dependent analysis when applied to many channels, 
is comparable with the precision of a conventional analysis which just uses the data from 
the muon pair asymmetry (there is some loss of statistical power due to the more restrictive 
requirements for events with accurate angular measurements). 

We have also demonstrated that the effect of IFI is adequately described by the leading 
order QED corrections, and that the asymmetry does vary greatly with the angular cut 
imposed, showing that, while the correction to a conventional analysis which integrates 
over all photon phase space is small, this is because of a large cancellation which requires 
respectful treatment. 

In experiments at proposed future electron-positron colliders [^TJ, the collisions be- 
tween the very dense bunches will produce radiation and lower the effective CM energy. 
This effect is similar to ISR, but depends not only on a standard QED radiator function 
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but also on the detailed bunch dynamics, which can vary from one collision to the next. 
This presents a serious challenge for conventional muon pair analysis at such machines, 
whereas this method is not disturbed by such a variation. 



9 Appendix 

The full expression for the probability density used in the unbinned likelihood fit has the 
following form: 



V( V ,z)^^e^\z\)e t (\z\)J2Az)p i ( S ,r ] )J2D i j (r ] )H j (z) (44) 

i j=l 

where Hj stand for different types of dependence upon z =cos6*: 



Hi = 1(1 + z 2 ) 

H 2 = ^(l + z 2 )llni±^ 

2 8 V ' Az \-z 

H?, Ha = z — In 

3 ' 2 \-z 

H 5 ,H 6 = z (45) 

H 7 , H 8 = z — ln - 

Az 1 — z 

H 9 = ^(l + z 2 )iln — 

9 8 V J 2 1- z 



Coefficients D\ -- D\ are expressed through constants A , A', B, defined in eqs. (p8|-p0|), 
and the ratios, \&i, ^2, ^3, of the convoluted functions f±: 

D[ = 1 

D\ = ¥M 

D\ = [A + A\s-M 2 z )/(2s)}B^\(r ] ) 

D\ = -A'B& 2 (ti) 

D\ = [A + A'(s-M 2 z )/(2s)] (46) 

D\ = -A'%( V ) 

D\ = [A + A'(s-M 2 z )/(2 s )]¥ 1 (r j ) 

Dl = -A'%( V ) 

Di = B¥ x (rj) 

where 



*<M = m, *M = l£M, %(n) = v 4M, (47) 

f+iv) f+(v) f+iv) 

with the horizontal bar denoting resolution smearing. The index i serves as a reminder 
that the resolution parameter a* used during the smearing, is different for different event 
classes i. Note that the functions ^\(r)) also implicitly depend on the acoplanarity cut £ - 
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The factor e^(|^|) stands for the selection probability of an event with a true acopla- 
narity £ true < £o when the cut is applied on the measured acoplanarity £ of the event. 
It was determined using the MCI sample, and is within ~ 10~ 3 of unity when a^ <C £o 
(which is the case for the barrel region \z\ < 0.7), decreasing slightly for larger z, where 
the £ resolution is worse. 

Similarly, the factor £t(|^|) takes into account the variation of efficiency with |cos#*| 
due to the cut on the missing transverse energy, eq. (p4|). It also was determined using 
the MCI sample, and decreases from ~ 96% in the barrel region to ~ 80% at the edge of 
the acceptance. 

e l (z) is the class-specific selection efficiency relative to the total efficiency for all classes, 
including those not used in present analysis, and is determined directly from the data 
(Fig. ^j. So are the functions p*(s, r/) which, for a particular class % at each initial energy 
point y/s, are essentially equal to the measured //-distributions, dominated by the j = 1 
term in eq. (f44"D , with small and calculable corrections from other terms in the sum. 

Finally, the normalisation constant M is determined from the condition that the total 
probability is equal to 1. 
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T data 
Z 141 MC 

n ttmc 




Figure 1: Distribution of selected class 1 muon pair events (as defined in subsection 5.2) 
vs missing transverse momentum (data points), compared to the Monte Carlo simulation 
of muon pair production (light histogram) and tau pair background (dark histogram). A 
cut xt < 0.1 accepts 96% of class 1 muon pair events and removes about 75 % of the 
background, which is further reduced by the tight acoplanarity cut. 
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Figure 2: Distribution of class |T| events from the central region |cos0*| < 0.2 with respect 
to the variables log(l/£) (a) and log(l/?y) (b). The positions of the peaks, shown by the 
arrows, determine the experimental resolutions in £ and r\. 
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Figure 3: Distribution of class p] events in the log(l/£) — log(l/?y) plane. Events from 
the lightly shaded triangular area above the diagonal have a high ISR probability, while 
the densely shaded area along the diagonal contains events with a high probability of 
significant FSR and/or IFI. The horizontal dashed line represents the cut on acoplanarity 
angle £o = 0.004, while the vertical dashed line corresponds to rj = 0.008. 
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Figure 4: cos ^'-distributions for the class 1 and class 2 events at Z peak, integrated over 
separate areas of the log(l/£) — log(l/?y) plane from fig. Q, together with fit results using 
the function (P5|): a) upper left quadrant, where the asymmetry is dominated by the linear 
term (d ~ 0); b) upper right quadrant, where the IFI-induced asymmetry is significant 
and positive (d > 0); c) lower left quadrant, where the IFI-induced asymmetry is large 
and negative {d < 0); d) the whole plane, where the IFI contribution is compatible with 
zero as a result of the cancellation. 
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Figure 5: Differential asymmetry defined in eq. (|26|) for the class 1 and class 2 events at 
Z peak, corresponding to four different areas of the log(l/£) — log(l/r/) plane from fig. [jf 
upper left quadrant (a), where the asymmetry is dominated by the linear term; upper right 
quadrant (b), where the positive IFI-induced asymmetry is dominant; lower left quadrant 
(c), where the IFI-induced asymmetry is large and negative; and the whole plane (d), 
with no significant asymmetry of any kind. Data points with error bars represent OPAL 
data, the solid histogram shows the MCI sample without the IFI contribution, while the 
dashed histogram shows the MC2 sample which contains IFI. 
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Figure 6: Typical dependence of the ratios \&i(^) = / '-(?]) / j +(77) (a), ^2(v) — 
Vf-iv)/ f+(v) (b) and ^3(1]) = i]f + (r])/ f+(rj) (c) on 77 before (dashed lines) and after 
(solid lines) convoluting each of the functions with the 77 resolution. In this example, the 
acoplanarity cut is £0 — 0.004, while the 77 resolution parameter is a„ = 0.002. 
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Figure 7: Selection efficiencies for classes 1 (dark histogram) and 2 (light histogram) 
as functions of cos 6', as determined from the sample of events passing all cuts for the 
asymmetry sample, including those on xt, £, V and cos 9*. 
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Figure 8: Fitted asymmetry in narrow slices of vs', with the initial-final radiation 
interference term fixed to its value measured in this analysis, B = 0.840. 
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